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Abstract 



We find a black hole solution with non-Abelian field in Brans-Dicke theory. It 
is an extension of non-Abelian black hole in general relativity. We discuss two 
non-Abelian fields: "SU(2)" Yang-Mills field with a mass (Proca field) and the 
SU(2)xSU(2) Skyrme field. In both cases, as in general relativity, there are 
two branches of solutions, i.e., two black hole solutions with the same horizon 
radius. Masses of both black holes are always smaller than those in general 
relativity. A cusp structure in the mass-horizon radius (Mg-rh) diagram, 
which is a typical symptom of stability change in catastrophe theory, does 
not appear in the Brans-Dicke frame but is found in the Einstein conformal 
frame. This suggests that catastrophe theory may be simply applied for a 
stability analysis as it is if we use the variables in the Einstein frame. We also 
discuss the effects of the Brans-Dicke scalar field on black hole structure. 
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I. INTRODUCTION 



For many years, there have been various efforts to find a theory of "everything" . The 
Kaluza-Klein theory was one of the candidates, which is constructed in a five dimensional 
spacetime. Jordan noticed in 1955 that in our four dimensional spacetime, a scalar field ap- 
pears by a compactification in the Kaluza-Klein theory and it gives a nonminimal coupling 
to gravity, meaning that this theory violates even the weak equivalence principle. Dicke 
thought that the weak equivalence principle must be guaranteed based on several experi- 
ments. Then, from the weak equivalence principle^ and Mach's Principle, which insists that 
an inertial force is determined by the distribution of matter all over the Universe, he and 
Brans constructed a new scalar-tensor theory, i.e., Brans-Dicke (BD) theory, in 1961 0. 
Since then the difference between general relativity (GR) and BD theory has been discussed 
in many aspects. Although BD theory itself is strongly constrained by several experiments 
(the BD parameter uj > 500), we believe that the theory may still be important from the 
following points of view: 

BD theory can be an effective field theory of a unified theory of fundamental forces. 
In particular the BD-type scalar field appears as a dilaton field in superstring theory. 

BD theory is one of the simplest extensions of GR. So if we wish to discuss something 
in a generalized theory of gravity, BD theory can be the best model to see a difference 
from GR. 

Moreover, a scalar field such as the BD scalar field may have an affect on many aspects 
in gravitational physics. For example, the inflationary scenario would be modified by an 
introduction of such a scalar field 0. Although the inflationary scenario was discussed 
originally in GR, since an appropriate model based on particle physics has not been found, 



^But for self-gravitating bodies, the weak equivalence principle is still violated (Nordtvedt effect) 
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it is important to recognize that the introduction of a scalar field can make a big change in 
scenarios of the very early universe. 

Black holes are also important in gravitational physics. We may expect that such a 
scalar field also affects some feathers of a black hole 0. However, since the gravity part in 
BD theory is conformally equivalent to that in GR, black hole solutions are not modified 
by the introduction of the BD scalar field for the case without matter or with a traceless 
matter field such as the electromagnetic field. As a result, for vacuum case or the case with 
the electromagnetic field, a conventional Kerr or Kerr-Newman black hole turns out to be 
a unique solution even in BD theory because of the black hole no-hair theorem in GR P|. 
Hence, here we shall discuss a non-Abelian black hole in BD theory, which has so far not 
been studied so much. For the case with the Yang-Mills field, however, we again find the 
same colored black hole as that in GR , because its energy-momentum tensor is traceless. 
Then, we discuss a "massive" non-Abelian field, i.e., a massive Yang-Mills (Proca) field, and 
the Skyrme field. We consider only the globally neutral case in this paper. 

After the introduction of basic ansatze and conditions in §.II, we present the Proca black 
hole solution and its properties in §.III. We find some difficulty in adopting catastrophe 
theory to the stability analysis. To resolve such a difficulty, we introduce new variables 
defined in the Einstein conformal frame in §.1V. We find quite similar properties of black 
hole solutions to those in GR: in particular a cusp structure appears in the mass-horizon 
radius diagram. This allows simple application of catastrophe theory in the stability analysis 
as it is. The effects of the BD scalar field on black hole structure are investigated in §.V. In 
§.V1, we discuss a Skyrme black hole, showing that its properties are quite similar to those 
in the Proca black hole. The concluding remarks will follow in §.V11. Throughout this paper 
we use units of c=h=l. Notations and definitions such as Christoffel symbols and curvatures 
follow Misner-Thorne- Wheeler [0]. 
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II. NON-ABELIAN BLACK HOLES IN BRANS-DICKE THEORY 



The action of BD theory is written as 



S = J d^x^ ^ (^(j)R - ^ V«0V"0j + 



(1) 



where = StcG with G being Newton's gravitational constant. The BD parameter is u, and 
Lm is the Lagrangian of the matter field. The dimensionless BD scalar field is normalized 
byG. 

For the BD field 0, the field equation becomes 

-2 



□c 



2uj + 3 ^ 



(2) 



Then, if the right hand side of this equation vanishes, that is, the energy-momentum tensor 
of the matter field is traceless, 0= constant turns out to be a solution, meaning that a black 
hole solution in GR is also a solution in BD theory. Hence, for the SU(2) Yang-Mills field, 
we find that the colored black hole 0] is a solution in BD theory too. Although we have no 
proof, we expect that for the case with a massless non-Abelian gauge field, no new type of 
black hole solution appears in BD theory. 

If a non-Abelian field is massive or effectively massive, however, = constant is no longer 
a solution. We will find a new type of black hole solution, and can discuss some differences 
from black hole solutions in GR. This is the reason for us to study a massive non-Abelian 
field here. 

We assume that a black hole is static and spherically symmetric, in which case the metric 
is written as 



1 



2Gm(r) 



e-^'^'^de + 



1 



2Gm{r) 



n -1 



(3) 



The boundary condition for a black hole solution at spatial infinity is0 

_ 2(2 + cj) 



lim m = M < oo, lim 6 = 0, lim 



3 + 2a; 



(4) 



^This choice of 0o at oo guarantees that G is Newton's gravitational constant B. 
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Note that a test particle far from a black hole does not move under an influence of this 
"mass" M, but feels a gravitational attractive force given by a gravitational mass Mg. Mg 
is defined from the asymptotic behavior of the time-time component of the metric and given 
as 

Mg = M + l:,\im {r8) . (5) 
For the existence of a regular event horizon, r/j, we have 

nih = m{rh) = 4 = S{rh) < oo . (6) 

We also require that no singularity exists outside the horizon, i.e., 

r 

m{r) < — for r > r^. (7) 

For our numerical calculation, we introduce the following dimensionless variables: 

r ^ r/rh, ffi ^ Gm/rh- (8) 

To write down the explicit equations of motion, we have to specify our models. In what 
follows, we discuss the Proca field and the Skyrme field, separately. 

III. PROCA BLACK HOLE 

We first consider a massive "SU(2)" Yang-Mills field (Proca field). The matter La- 
grangian is now 



TrF' - -^TrA' , (9) 



where gc and /i are the gauge coupling constant and the mass of the Proca field, respectively. 
F is the field strength expressed by its potential A as F — dA -\-A/\A. For the spherically 

symmetric case, we can write the vector potential as 

hir t) 

A — air, t)Trdt H ^ — Trdr + {d{r, t)T0 — [1 + w{r, t)\Ts}d9 

r 

+{[1 + w{r, t)]Te + d{r, t)T^} sin Odcj) (10) 



as Witten showed 0, where Tr,T0, and r,/, are the generators of su(2) Lie algebra. We 
adopt the 't Hooft ansatz, i.e., a = 0, which means that only a magnetic component of the 
Proca field exists. We also set b = 0.^ In the static case, we can set d = 0. Now, our 
potential is 



A = [l+w{r)]{-T^d9 + Tesm 



(11) 



The boundary condition of the Proca field for its total energy to be finite is 



lim w = —I 

r— >oo 



(12) 



We define dimensionless parameters as 



11 = fi/ {gcmp) , Xh = rh/ {Ipl Qc) 



(13) 



Ip = G^/^ and rup = G^^/^ are the Planck length and mass defined by Newton's gravitational 
constant, respectively. 

Under the above ansatze, we find the following basic equations: 



dm / , d(j) 
dr \ dr 



2mx 

f / 1 26 \ df 



^ 2 fdwV 



Xl \ dr I \ dr 



+ 



1 



X2 



d) dr 



2cu + 3 {2uj + 3)0 df 



(14) 



d5 r / ^ ^ dcp' 
df 6 \ df 



fd4>Y 4 fdw\ 



-1 



2^ + .-^ (1 



2m 



r 



1 r d(f) 
Xl^df 



1_^2\2 ,,,, , 1^,2X2 



4(cj + l)fx^Xi fl+w 
2cu + 3 



2 / m\ d(f) 4/2^ /I + w'^ ^' 
f \ f J df 2uj + 3 \ f 



(15) 



d 



dr^ (p \ dr I \ r 



1 r d(j) 
Xl4)df 



l-W^V . 4(c^ + l)/i2A2 + 



+ 



2cu + 3 



^If the Yang-Mills field is massless, we can always impose this condition via a gauge transformation. 
In the present case, however, we can just show that this is consistent with the field equation. 
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1-- 



4^2 



r J dr 



d^w 1 dw d(j) / 2m 
df"^ d) df df \ f 



1 r dw 
Xl^df 



2uj + 3 
1 



1 + w 



(16) 



:i+w){fi'xi 



w{l — w) 



2m dw 

f2 




+ 



4(cj + 



(17) 



{2uj + 3) 

As for the boundary condition at the event horizon, in order for the horizon to be regular, 
the square brackets in ([I5|), (|T6|) and (|1^) must vanish at the horizon f = 1. Hence, we find 



that 



dw 
df 

df 



f=i 



0feA|(2^ + 3)(1 + Wh){[i^\l - Wh{\ - wh)} 
(2u; + 3){(1 - wl)-^ - 0,A|} + + 1)A2(1 + w^Y ' 

4/i2(l + U;;,)2A|0;, 



f=l 



(18) 
(19) 



(2^ + 3){(1 - «;2)2 _ + 4/i2(c. + l)A^(l + w^Y ' 

where Wh = w{rh) and (ph = (t>{rh). As a result, if/i and 0/i turn out to be shooting 
parameters and should be determined iteratively so that the boundary conditions @) and 
(P^) are satisfied. 

In Fig.l, we present a numerical solution with = 0.5lp/gc and /i = 0.15gcmp^ We 
set uj = 0. Although this is not consistent with the present limit from experiments, we 
choose this value because we wish to clarify the difference from black holes in GR. For a 
massive non-Abelian field, the node number of the potential w{r) is limited by some finite 
integer. Here the node number is chosen to be the smallest value, i.e., one. The dotted line 
denotes the Proca black hole in GR ^ with the same parameters, i.e. Vh = 0.5lp/gc and 
/i = O.lbgcmp, which we show reference. 



^As we see later, for fixing r^, we have two solutions in BD theory as in GR. Here, we show the 
field distributions in the solid-line branch in Fig.|3|. 
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FIG. 1. The solution of the Proca black hole in BD theory with lo = for = 0.5lp/gc 
and /i = 0.15gcmp{ (a) (j){r), (b) fh{r), (c) w(r), (d) 6{r) ). : The Proca black hole in GR is also 
depicted as a reference by a dotted line. The arrow in (c) shows the Compton wavelength of the 
Proca field (1/^). 

As seen from Fig. 1(a), the BD scalar field decreases monotonically as 

1 + -) , (20) 



where Mg is a constant and called the scalar mass |]T0[. Fig. 1(c) shows that the non-trivial 



structure of the non-Abelian field extends to the scale of the Compton wavelength of the 
Proca field (~ 1/yu), which is shown by an arrow. From Fig. 1(d), you may not see a clear 
difference between the lapse function 6 in BD theory and that in GR, but 6 in BD theory 
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falls as 1/r as r — s> oo while S in GR vanishes much faster than 1/r. In fact, from Eq. ([I5D 



we find 



d6 1 d(f) 



2GM, 



s 



(21) 



dr (j) dr 



near spatial infinity. This gives the relation between M and Mg as 



9 ~ 



M + 2M 



s 



(22) 



To see a property of a family of black hole solutions, we show the relation between the 
gravitational mass Mg and horizon radius r/j in Fig.2. The dotted lines denote the Proca 
black hole in GR with the same parameters, i.e. fi = O.lgcnip or fi = O.lbgcrrip, and the 
dot-dashed lines represent the Schwarzschild and colored black holes, respectively, which we 
show as references. 

As we mentioned in Fig. 1(c), the non-trivial structure of the non-Abelian field is as large 
as the scale of the Compton wavelength (~ l//i)- This is responsible for the existence of 
a maximum horizon radius (~ l//i) as in GR. That is, beyond this critical horizon radius, 
a non-trivial structure is swallowed into the horizon and then cannot exist, resulting in a 
trivial Schwarzschild spacetime. 

The Schwarzschild black hole is a trivial solution (m = Mg, 6 = 0, (p = (po and w = — 1), 
which has no upper bound for a mass or a horizon radius. If the YM field is massless, a 
family of colored black holes also exists as a non-trivial black hole, where the BD scalar field 
is = 00 = constant. There is also no upper bound for horizon radius as in GR. 



9 



FIG. 2. Mg — rfi diagram of the Proca black holes. The solid loop lines denote Proca black 
holes with /j, = O.lgcmp and OAbgcTUp in BD theory (to = 0). We depict those in GR with the 
same parameters by dotted lines. The Schwarzschild and colored black holes are also shown as 
references. 

The mass of the Proca black hole in BD theory is always smaller than that in GR (see 
also Fig. 1(b)). This is just because the value of the BD scalar field near the black hole 
is larger than that at infinity, which means the effective gravitational constant is always 
smaller than G. Therefore the mass concentration by gravitational attractive force may get 
smaller. 

In GR, there are two branches of black hole solutions: One is stable and the other is 
unstable. Those two branches coincide at a critical horizon radius or at a critical mass, 
where we find a cusp structure on the gravitational mass Mg- horizon radius Vh diagram. 
This cusp structure is a typical symptom of stability change in catastrophe theory ||1T|. The 
stabiUty analysis by catastrophe theory agrees with that by linear perturbations. 
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FIG. 3. Mg-rfi diagram for Proca black holes in BD theory [lo = 0) and in GR. The mass of 
the Proca field is /U = Q.lbgcmp. 



To see the detail and compare our solution with that in GR, we depict the enlarged 
diagram for /i = 0.15(7c'^p in No cusp structure appears in BD theory. Though the 

solution curves seems to merge at the point D, we find two solutions exist at D, which can 
be distinguished from field distributions. In GR, the maximum points of horizon radius Th 
and of gravitational mass Mg are the same, i.e., the point C in Fig.|^. In BD theory, however 
those two points, A (maximum horizon radius) and B (maximum gravitational mass), are 
different. This result does not depend on the choice of the BD parameter uo and the mass 
of the Proca field /i. In particular, a cusp structure disappears in BD theory as mentioned 
above. One may wonder whether catastrophe theory can be simply applied to stability 
analysis as it is in BD theory, although it works quite well in GR |Il2|. For fixing r/j, we 
still have two solutions in BD theory as in GR. Is there any correspondence of those two 
solutions to two branches in GR? We expect that there are two types of black holes in the 
BD theory as well. 



As we discussed in our previous papers |]T2|, if we divide the total energy density ptotai 
into a kinetic term pp2 and a mass term p^2, one of the main differences between the 
two branches in GR (the solid-line and the dotted-line branches^ in Fig.^) comes from the 
difference of dominant ingredient, i.e., in the solid-line branch, is dominant compared 



^In 1 12], we called them high-entropy and low-entropy branches, respectively. 
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D 



FIG. 4. Distributions of the energy density for (a) the soHd-Hne and (b) the dotted-hne branches 
of the Proca black holes with /i = O.lbgcmp both in BD theory (uj = 0) and in GR. The horizon 
radii of the black holes are = O.Ol/p/^c- 



In the solid-line branch, the black hole and non-Abelian structure are rather independent. 
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D 



D 

DD 



FIG. 5. rfi-Mg (scalar mass) diagram of the Proca black holes with fi = O.lS^cW-p in BD theory 
(w = 0). The dotted and solid lines correspond to those in Fig.^. 
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(b) 

FIG. 6. (a) Mg-l/T diagram for Proca black holes with /x = OAbgcirip in BD theory (a; = 0) 
and in GR. 

FIG. 7. (b) rfi-Bfi diagram for Proca black holes with /x = O.lbgcmp in BD theory (a; = 0) and 
in GR. 

IV. PROCA BLACK HOLE IN THE EINSTEIN CONFORMAL FRAME 

The gravity part of BD theory is conformally equivalent to that of GR, and a description 
by use of the Einstein conformal frame sometimes gives us simpler basic equations and easier 
analysis because the coupling of the BD scalar to gravity is moved to a matter term and 
the gravity part is just described as in the Einstein frame, which is already familiar. Hence, 
here we shall reanalyze our present problem in the Einstein conformal frame. We consider 
a conformal transformation 

(f) 

gab = -rQab ■ (29) 

00 

The equivalent action S = S/(f)o is given as 

^TrF^ + — 



5* = y d'^x^J^-§ 



-R - - Va(^V> - — — ^TrF' + exp(-/^/3(^)TrA 



2 



Kp" I 00. 



, (30) 



.-^mm, ,^C-^r'. (31) 
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For a black hole solution, if we define spherically symmetric coordinates in the Einstein 
frame as 



00 



1 



2Gm(r) 



2Gm{r) 



n -1 



rff2 + f'^dn^ 



(32) 



we find 



M = lim m(f) = Mq — Ms, fh = rh\ 



(33) 



where variables with a caret ' denote those in the Einstein frame. We also introduce dimen- 
sionless variables and parameters as 



r = r/r-h , m = Gm/rh , K = rhQc/lp , ^ = ^/ 



m 



p ■ 



(34) 



The basic equations are now 

2 



dm 
df 



d5_ 

df 

d^ 
df^ 



d'^w 
df'^ 



1 



h 

+ ^exp(-v/ 
— 47rf ^ 



df I \ df , 



^ 2fh 

f 



20oA 



^7r/5(^)(l + w) 



dw 



df 



+ 1 



2m 

f 



(j)oXlf V dr 

^ ' 'dip' 
df 



(35) 
(36) 



2m 
r 



-— + 



1 — w 



2\ 2 



0oA^ 



+ 2f§^ exp(- 78^/3 v?) ^^^"^ 



^0 



r 



1 



47r0o 
2m 
r J 

w{l — w"^) 
r 



exp 
1 



y8^/3<^) 
diu j 2m 

+ /i^A|exp(-A 



1 + 



r 
f 



(37) 



/8^/3(^)(l + w) 



+ 2r — exp 

00 



l+W 



(St 



As we expected, these are simpler than those described in the BD frame. The boundary 
conditions are similar to the ones in the BD frame. From our numerical calculation, we can 
show that M = Mg because S vanishes faster than f~^. 

First, in Fig.|^, we show the Mg-fh diagram in the Einstein frame that is related to Fig.|^ 
by conformal transformation. Surprisingly, we recover a cusp structure even in BD theory. 
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D 



(b) 

FIG. 8. (a) Mg-rfi diagram in the Einstein frame for Proca black holes in BD theory (cj = 0) 
and in GR. The mass of the Proca field is ;U = O.lBgcmp. We find a cusp structure, which indicates 
a stability change via catastrophe theory. 

FIG. 9. (b) Mg-l/T diagram in the Einstein frame for Proca black holes in BD theory (uj = 0) 
and in GR. The mass of the Proca field is /.i = OAbgcTrip. 

This suggests that catastrophe theory will be simply applied in a stability analysis for 
non-Abelian black holes in BD theory as well. 

From the point of view of catastrophe theory stability changes at a cusp point in the 
control parameter-potential function diagram. In GR, if we regard gravitational mass and 
black hole entropy (or equivalently the area of the event horizon) as a control parameter and 
a potential function, respectively, we find a cusp C in the Mg-r^ {Mg-Vh) diagram (Figs.|,||), 
which is a symptom of stability change in catastrophe theory. In fact the stability of the 
black hole does change at this cusp point C. In BD theory, however, a cusp structure does 
not appear in Fig.^, while it does in Fig.|[ This suggests that if we use the variables in 
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the Einstein frame, we can simply apply catastrophe theory to the stability analysis in BD 
theory as it is. From Fig.^, catastrophe theory predicts that stability change can occur at 
the point A. From Fig.|], however, no such prediction is possible. 

To study stability, we have another method, i.e., a turning point method for thermody- 
namical variables |13|. Stability will change at the point where d{l/T)/dM = oo. In GR, 



we understand that a stability change occurs at the point C in Figs.^p. This is consistent 
with analysis by catastrophe theory. In BD theory, d{l/T)/dM = oo occurs at the point B 
in Fig.|^ (BD frame), which is inconsistent with the stability analysis by catastrophe theory. 
However, if we use the Mg-l/T diagram in the Einstein frame, the divergence occurs at 
the point A, which is consistent with catastrophe theory. To understand this inconsistency, 
we have to remember that variables in the turning point method should satisfy thermody- 
namical laws, in particular the "mass" of a black hole should satisfy the first law of black 
hole thermodynamics. In fact, the gravitational mass in the BD frame does not satisfy the 
first law of black hole thermodynamics, while it does so for the variables in the Einstein 
frame (Fig.P)|, and therefore the turning point method could be applied. We expect that a 
stability change occurs at the point A but not at the point B, and this is consistent with 
catastrophe theory. 

These conjectures for stability should be justified by analyzing linear perturbations of 
black holes and black hole thermodynamics 0]. 

V. THE EFFECTS OF THE BRANS-DICKE SCALAR FIELD 

By use of the Einstein frame, we also understand easily a qualitative difference between 
black holes in BD theory and in GR. As we see in the action (|30D , the coupling of the 
BD scalar field appears in the mass term. This coupling reduces effectively the mass of 
the Proca field by a factor exp{—K,f3(f/2) because ip is monotonically decreasing to zero as 
r — i> oo (see Fig.|l](a)). In GR, as the mass is reduced, the Proca black hole shifts in the 



^We can show that thermodynamical variables in the Einstein frame satisfy the first law. |14| 
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FIG. 10. (a) Mg-rh diagram in the BD frame (b) Mg-fh diagram in the Einstein frame for 
several values of uj. The mass of the Proca field is ^ = 0.15gcmp. 
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FIG. 12. (a) Mg-rh diagram and (b) Adg-l/T diagram normalized by g'^ = ^/(po gc in the 
Einstein frame for several values of uj. The mass of the Proca field is /i = O.lbgcmp. 

From Fig.|l^, in which the effect is absorbed in normahzation by g'^, we find the deviation 
from GR is quite similar to the behavior when changing a mass of the Proca field in GR 
(see Fig.0). This means that effects (1) and (2) really explain the deviation from GR. 

In Fig.|l^, we depict the gravitational mass and the lapse function in terms of u for 
fixed horizon radius {vh = 0.5lp/gc) and fixed mass of the Proca field {fi = O.lbgcTrip). The 
solid and dotted lines correspond to those in the solid-line and in the dotted-line branches, 
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FIG. 13. oj dependence on the gravitational mass Mg and the lapse function 5 of the Proca 
black hole in the BD frame for a fixed = 0.5lp/gc and = O.lbgcmp ((a) , (c) : uj < 1 and (b) 
, (d) : a; > 1). The dotted and solid lines correspond to the dotted- and solid-line branches in 
Figjl. 
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FIG. 14. oj dependence on the gravitational mass Mg and the lapse function 5 of the Proca 
black hole in the Einstein frame for a fixed = 0.5lp/gc and /i = 0.15gcmp ((a) , (c) : a; < 1 and 
(b) , (d) : uj > 1). The dotted and solid lines correspond to the dotted- and solid-line branches in 
Figjl. 

We can see that the gravitational mass approaches some constants as oo, which 

correspond to those in GR. In the Einstein frame, the horizon radii in both branches approach 
the Schwarzschild radius in the limit of u; = —3/2, resulting in a trivial Schwarzschild black 
hole(Fig.|T|). 

This is because the matter contribution will vanish in this limit as we discussed above 
(00^ 0). In the BD frame, however, we find that the dotted-line branch changes faster 
than the solid-line branch and both horizon radii in the limit of u = —3/2 are different from 
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the Schwarzschild radius. Then non-trivial black holes can exist even for uo = —3/2. This is 
consistent with the above result in the Einstein frame because the conformal transformation 
becomes singular for u = —3/2. 

Although any value of u; > —3/2 does not give a ghost, we find a negative mass con- 
tribution in the BD frame, resulting in that M becomes negative for u < {ucr < — 1). We 



show m(r) for several values of u in Fig.[T5|. This does not mean, however, that we have a 
negative-mass black hole, because the gravitational mass Mg is still positive. A test particle 
moving around a black hole feels an attractive force given by Mg, which is always positive. 
The effect of negative M could be observed in a time delay, which changes its sign for 



negative M |]T0 



In the Einstein frame, m(r) is monotonically increasing as r ^ oo, resulting in a positive 
mass M, which is the same as the gravitational mass Mg (Fig.0). As we know, in BD 



theory, we can define several masses |jTO|] . The reason why we have several masses is because 
the BD scalar field decreases as r~^, which is responsible for having different masses in each 
frame, and the scalar field itself also gives a contribution into a mass energy as a scalar 
mass Mg. In the vacuum case, we find a negative M for u < —1. In our case, however, the 
BD scalar field is concentrated by the gravitational attractive force of the black hole. This 
changes the sensitivity s just as for a self-gravitating star. From the asymptotic behavior of 
Qoo and grr, we find a relation between M and Mg as 



Mg _ 00 + 2- s 



(39) 



M uj + l + s 

where s is a sensitivity (see Eq.(11.83) in [|lOl)- Then the sensitivity s could be evaluated as 

{2uj + 3){Mg - M) 



s 



(40) 



2 2(M + Mg) 

For a Schwarzschild black hole (M = Mg), s = 1/2. If u = —3/2, however, s = 1/2 even 



if M 7^ Mg. We show the sensitivity s in Fig.p!7|. From Eq.(^), M becomes negative for 
uj < — (1 + s) (< —1). Then for a given uj (< —1), M of the Proca black hole with smaller 
sensitivity than Scr = —{1+uj) becomes negative. It may correspond to smaller black holes in 
the solid-line branch from Fig.^. When lj oo, we find M = Mg but s ^ 1/2. The reason 
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□ 



□ 



□ 



FIG. 17. rh-s diagram for several values of uj. The mass of the Proca field is /i = ClS^c^p 



The dotted and solid lines correspond to the dotted- and solid-line branches in Fig.lC 



23 



FIG. 18. Lu dependence on the sensitivity s of the Proca black hole for a fixed rh = 0-5lp/gc 
and fi = O.lBgcmp {{a) u> < 1 and (b) a; > 1). The dotted and solid lines correspond to the dotted- 
and solid-line branches in Fig.p!o[ 



VI. SKYRME BLACK HOLE 

In GR, non-trivial black holes with a massive non-Abelian field have quite similar prop- 
erties, which we classified as Type II in [|T2|. How about black holes in BD theory? To see 
whether the above results for the Proca black hole are generic or not, we shall study the 
Skyrme field as another example of a "massive" non-Abelian field. 



The action of the Skyrme field Lm is SU (2) x SU (2) invariant and given as |jT5 



2 



TrF^ - ^Tr , (41) 



?>2gl— 4 

where fs and Qs are coupling constants. Qs is related to Qc for the Proca field as 



Qs = V4ngc . (42) 

The "mass" parameter of the Skyrme field, /i, is defined by yU = fsOs- F and A are the field 
strength and its potential, respectively. They are described by the 5'f7(2)-valued function 
U as 

F = A A A, A = U^WU . (43) 
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In the spherically symmetric and static case, we can set U, 



U{x) = cosxi^) + i sin x{r)cr if ^ 



(44) 



where cTj and are the Pauli spin matrices and a radial normal, respectively. The boundary 
condition for the total field energy to be finite is 



lim X = . 



(45) 



For simplicity, we solve the present system in the Einstein frame. The equivalent action 

S = S/cpo is 



S = J d'^x\J^-^ 



With the dimensionless parameter fg = fsl'^-p-, the basic equations are now 



dm 
df 



^ = 2n 



? 1 - Hi 

\ r 



df / 



'dxVfl 
df \ dn 



exp 



sin^ X 



+ sin^ X 2—^ exp 



47r0oA^r' 



d5 
df 

d^ 
df^ 



# f sin^ X fs . /^n -^ , sin^ X 
Stt— exp (- V87r/3(^) + —j—2 



= 1- 



2m' 



(if [ f 
oirp— exp (— V 87rp(/)) 



2m 



f df 200 \ '^'^ , 

exp (-v/87r/3(^) / 

Anf'^Xlf^ exp (-A/87r/5(^) + 2 sin^ x \ 
^ '" sin2x 



2^2 r2d^ c?x 



(46) 



(47) 
(48) 



(49) 
dx 



^-^Kfs^ - 47rV87r/3r A;,/, + sm2x — 

dr dr ' 



/ 2m 

+ 1-^ 
\ r 



47rf A|/2 exp {-VS^pip) + 2 sin^ x 



47rA2/2exp(-v^/?^) + 



df 

sin^ X ^ 
f^ , 



+ 



(if f 



8vr^ exp (-78^/3^) + ^^ 



2m 



(50) 



As in the case of the Proca black hole, the square brackets in ( ^91 ) and (|50| ) must vanish at 
r/j for the horizon to be regular. Hence 
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D □ 
□ 



FIG. 19. The solution of the Skyrme black hole with fh = l.Olp/gc and fs = O.OSmp in BD 
theory ( (a) (j){r), (b) m{r), (c) x(^)) (d) ^{f) ). : The Skyrme black hole in GR is also depicted 
as a reference. 

Here we set the parameters as 
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fh = l.Olp/gc, fs = O.OSmp, u = -1.4 . (53) 
The dotted lines are those in GR with = l.Olp/gc and fs = O.OSmp. The solutions 



correspond to solid-line in Fig.|20|. We have shown only for a solution with one node number. 
For a Skyrme black hole, rather than the node number, the solution is characterized by the 
"winding" number defined by[] 

Wn = -\xh - x{oo) - sm{xh)\ ■ (54) 

We show for a solution with the "winding" number one. Note that the comparison is 
made in the Einstein frame for a fixed fh, which does not mean the horizon radii with 
different u in the BD frame are the same. 

S falls faster than because (^8]) is 



as f ^ cxD and (f vanishes faster than (see Figs.[l9|(a), (d)). Then, as in the Proca black 
hole, M = Mg. 

To study the properties of a family of black holes, we depict the Mg-fh and Mg-Vh (the 
BD frame) diagrams in Fig.|20| and the Mg-l/T and Mg-l/T diagrams in Fig.pT. 



^ For a particlelike solution (Skyrmion), the value of x the origin must be vrn, where n is an 
integer and \n\ denotes the winding number of the Skyrmion. In the case of the black hole solution, 



it is topologically trivial. But Wn defined by (54) is close to n, so we shall also call it the "winding" 
number. 
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FIG. 21. (a) Mg-l/T diagram in the BD frame and (b) Mg-l/T diagram in the Einstein frame 
for Skyrme black holes. We set fg = O.OSmp. 

We find that the results are quite similar to those for the Proca black holes. We have 
a cusp structure in the Mg-fh diagram in the Einstein frame, but it disappears in the BD 
frame. 

Most properties found for the Proca black hole apply to the Skyrme black hole as well. 
This suggests that a universal picture for non-trivial black holes with massive non-Abelian 
fields is possible. 
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VII. CONCLUDING REMARKS 



We have analyzed non-Abelian black holes (Proca and Skyrme black holes) in BD theory 
and shown some differences from those in GR. The Einstein conformal frame makes our 
analysis easier. The effect of the BD scalar field can be reduced into two parts in the Einstein 
frame: the effective change of mass of the non-Abelian field, i.e. /i —>■ fiexp{—Kl3ip/2) or 
fs fsQ^Y>{.~i^Pvl'^)i s-iid the renormalized coupling Qc \Rk)gc oi Qs \/(hgs and 
fs ~^ fs/V^Po- As a result, the solutions shift in the left-upper direction in the Mg-r^ 
diagram. Although we recover the Schwarzschild black hole in the limit of ^ — | in the 
Einstein frame, we still have non-trivial black holes in the BD frame in the same limit, 
because the conformal transformation becomes singular then. 

Secondly, we have analyzed for various values of u. When u > 500, the difference 
from GR is so small that we will not see any observational difference. The solutions for 
^1 ^ — 1 seem to be somewhat pathological, because the mass function m becomes 
negative in the BD frame, resulting in negative value of M. However, even in such cases, Mg 
is always positive, therefore a test particle around such a black hole still feels an attractive 
force. 

Thirdly, we find that the cusp structure in Mg-Vh diagram does not appear in the BD 
theory although it was found in GR and provided us a new method for stability analysis 
via catastrophe theory, while it exists in the Einstein frame. This suggests that a stability 
change occurs at a cusp point in the Einstein frame. The justification of this conjecture and 
the proper analysis including that by linear perturbations will be given elsewhere [0]. 

In this paper, we have studied a globally neutral type of non-Abelian black holes in 
BD theory. A globally charged black hole, i.e., a monopole black hole may be much more 
interesting. That is because charged black holes are important in the context of cosmology, 
in particular, in the relation with a dynamical monopole (topological inflation) [|l^ . This is 
under investigation. 
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